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In order to study the interplay between Kondo and Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction, we calculate the spin-spin correlation functions between two Kondo impurities coupled 
to different sites of a half-filled open Hubbard chain. Using the density-matrix renormalization 
group (DMRG), we re-examine the exponents for the power-law decay of the correlation function 
between the two impurity spins as a function of the antiferromagnetic coupling J, the Hubbard 
interaction U and the distance 7? between the impurities. The exponents for finite systems obtained 
in this work deviate from previously published DMRG calculations. We furthermore show that 
the long-distance behavior of the exponents is the same for impurities coupled to the bulk or to 
both ends of the chain. We note that a universal exponent for the asymptotic behavior cannot be 
extracted from these finite-size systems with open boundary conditions. 



I. INTRODUCTION 

The Kondo effeciji is one of the oldest and most stud- 
ied correlation phenomena in condensed matter physics. 
It has regained vital interest by single-impurity setups 
and microscopic measurements i^"— One of the present 
main foci is the extension of the Kondo cloud and its 
experimental measurement i^^— The basic idea is that the 
conduction band electrons will form a highly correlated 
quantum state with the impurity spin and screen iti^ This 
results in nonzero spin-spin correlations between the im- 
purity spin and the conduction band electrons. The de- 
cay of these spin-spin correlations has been the subject 
of many theoretical studiesjiSrJ^ In particular, the corre- 
sponding envelope of these correlations was found to cross 
over from a 1 /R decay to a 1 /B? decay at the Kondo co- 
herence length, where R denotes the distance between 
the impurity and the conduction band siteJ^ 

If a second impurity is added to the system, one 
will have two competing interactions: the Kondo ef- 
fect and the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction^^"— The Kondo effect will screen the im- 
purity spins individually, leading to vanishing spin-spin 
correlations between the two impurity spins. The RKKY 
interaction favors a magnetic interaction between the im- 
purity spins, i.e., one observes strong correlations be- 
tween the two impurity spins. The interplay between 
these effects in the two-impurity problem has attracted 
much attentionii2ri2^ 

In this work, we want to focus on the static spin-spin 
correlations between two Kondo impurities. Another as- 
pect, which is usually neglected, is the presence of cor- 
relations in the conduction band system. They can in 
principle further modify exponents and also introduce 
additional functional dependencies, such as, for exam- 
ple, logarithmic corrections. Therefore, we examine the 
system in the presence of a finite U in the conduction 
chain. As an analytical treatment of the two-impurity 



Kondo problem is not available, one has to rely on numer- 
ical solutions. We will use the density-matrix renormal- 
ization group>2ii^ (DMRG) to obtain the ground state 
and to calculate the spin-spin correlation functions. In 
particular, we are interested in the form and decay of 
these correlations. Hallberg and Egger— calculated ex- 
actly these spin-spin correlations shortly after the de- 
velopment of the density-matrix renormalization group. 
From their data, they have argued that the correlations 
will show a power-law behavior in the long-distance limit 
(S/S//)i^ (X \/R^ for two Kondo impurities irrespective 
of the interaction U in the chain. 



Based on refined numerics for larger systems, which are 
accessible due to the great increase in computer power, 
we re-examine these exponents. For [/ > 0, we show that 
even with this increase in computational resources, one 
cannot easily identify a simple 1/R^ power-law behavior 
for the spin-spin correlations between two Kondo impu- 
rities. Only if we include a J7-dependent logarithmic cor- 
rection, our data are compatible with an exponent a ~ 2. 
This investigation shows that even for system sizes acces- 
sible today, an unbiased estimation of exponents for these 
long-range correlations is very difficult, especially if one 
has to expect logarithmic corrections. It is then impor- 
tant to understand how well such numerical calculations 
can reveal exponents expected to rule the decay of cor- 
relation functions for sufficiently large distances of the 
impurities. 



The paper is organized as follows. After presenting 
the model and a brief discussion of the method and its 
problems in Sec. |TT1 we present our results in Sec. IIIIl 
starting with impurities attached to the chain ends. The 
main results and conclusions of the paper are summarized 
in Sec. HVl 



II. MODEL AND METHOD 

We study two spin-1/2 Kondo impurities attached to 
a one-dimensional Hubbard chain. Figure[I] shows the 
corresponding setup. The Hamiltonian 

H = Hc + Had 

can be divided into two parts. The first contribution 

Z,CT i 

models the Hubbard conduction band with a Coulomb 
repulsion [/ > at each site. Here denotes the usual 
fcrmionic annihilation (creation) operator at site i and 
rii represents the particle number operator. Through- 
out our work, we adopt open boundary conditions and 
the hopping parameter is set to t ~ 1. For [/ = the 
conduction band He reduces to a tight-binding model. 

The coupling of the two Kondo spins to the conduction 
electrons is modeled by the s-d exchange term 

Had = — .//S/ • s„ — JllSjI ■ S,n, 

where S/i"//) {S'^,S'^,S^) is the spin operator 

of the Kondo impurity which is attached to site 
n(m) of the Hubbard chain. The operator s„(,„) = 
(*n(m)'*^(m)'*n(m)) the Spin Operator at the conduc- 
tion band site. We restrict our study to the case of 
antiferromagnetic coupling constants of equal strength 
Jj = Jjj = J < 0. 

We use a standard density-matrix renormalization 
group (DMRG) algorithm^"— for open boundary con- 
ditions to carry out the calculations. The setup with the 
two impurities attached to the ends of the chain is compu- 
tationally less demanding and hence more accurate than 
the analysis of the bulk limit [R <^ -^/2). Since it will 
turn out that the results in the long-distance or strong- 
coupling limit are the same for these two setups, the ma- 
jority of our calculations is performed with the impuri- 
ties coupled to the ends. This is still very costly because 
for every different inter-impurity distance, a new DMRG 
setup as well as a new DMRG ground-state calculation 
is needed. Furthermore, extracting the true long-range 
behavior of correlation functions within the DMRG is 
known to be difficult since all correlations are either long 
ranged or purely exponentially decaying, and even for 
power-law correlations this modeling only works for not 
too long distances (cf. Rcf. onlincciteSchollwock2005). 
In order to estimate the minimal number of basis states m 
that has to be kept within the DMRG, the upper panel of 
Fig. [2] shows the end-to-end correlation function (sfsf^) 
for a half-filled Hubbard chain with U = 1 and differ- 
ent values of m. At first sight, the DMRG data for the 
correlation function itself looks the same irrespective of 
the number of kept states. However, a more thorough 
analysis by means of double-logarithmic central differ- 
ences (to be introduced later, see Eq. ([5])), shown in the 
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FIG. 1. Two spm-1/2 impurities (filled circles), separated by 
the inter-impurity distance R — \m ~ n\, are attached to a 
one-dimensional Hubbard chain of length L. 
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FIG. 2. (Color online) Upper panel: Comparison of DMRG 
data for the end-to-end correlation function (sfsf,) of a half- 
filled Hubbard chain with U = 1 on a linear scale for different 
numbers m of kept states. Lower panel: Double-logarithmic 
central differences giving the exponent a{R) extracted from 
(sfsf,) for f/ = 1 in the spirit of Eq. ((4]). 



lower panel of Fig. ^ reveals that one needs to choose a 
very large number of basis states m in order to obtain 
a reasonably converged estimate of the long-distance de- 
cay respectively the associated power law. Thus, our first 
conclusion is that the question of convergence very subtly 
depends on the property one is eventually interested in. 
From Fig. [2] we also conclude that even a reliable extrac- 
tion of the exponent is very difficult for longer distances. 
Consequently, in our DMRG calculations, we typically 
keep m = 1300 basis states for chains up to lengths of 
L ^ 120. 

At half-filling, we learned from exact diagonalization 
results for small systems of even length L that the ground 
state is nondegenerate for S^^f. = 0, if we attach one 
Kondo impurity to an even site and the other one to 
an odd site. The inter-impurity distance is given by 
R ~ |m — n|. The DMRG calculations are carried out 
in the two Abelian U(l) symmetry sectors defined by the 
electron number TV and the z-component of the total spin 
S^^t = Sf + SI J + S^. Introducing full SU(2) symmetry 
in the spin sector can reduce the computational effort 
slightly, but will not alter our main observations. 
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III. CORRELATION FUNCTIONS 

The main focus of our work is the analysis of the spatial 
behavior of the spin-spin correlation functions (S/S//)^, 
between the two attached Kondo spins at half-filling and 
zero temperature. In particular, we want to determine 
whether for infinitely large distances R a power law of 
the form 

(S/Sj,)fl cx i?-" (1) 
with a constant exponent a ^ a(R) exists. 

A. Impurities coupled to the ends 
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The impurities are attached to the first and the last 
sites of a Hubbard chain of length L. The distance be- 
tween the impurities therefore is i? = L — 1. The corre- 
lation function is 2kF oscillatory, with kp given by the 
filling, e.g., kp = 7r/2 for half-filling. We will concentrate 
on odd distances R to avoid these oscillations. 



1. Non-interacting conduction chain U = 

We start by re-examining the system with a nonin- 
teracting conduction chain. For \J\ oo, the Kondo 
impurities form a rigidly bound singlet with the spins 
at the conduction band sites they are attached to, i.e., 
the first and the last sites. In this case, the system de- 
couples into three parts (cf. Ref. [H): the two singlets 
and an effective chain of length L — 2 in-between. Thus, 
for strong couplings |J| ^ 1, there is almost no hop- 
ping between sites 1 and 2 as well as between L — 1 and 
L. These considerations form the basis of a perturbative 
treatment which has originally been proposed for a single 
Kondo impurity coupled to one end of a noninteracting 
conduction bandi^^ In Appendix [Xl we extend this result 
to the case of two impurities attached to the ends. Our 
DMRG results can therefore be directly compared with 
perturbation theory for strong couplings |J| ;3> 1 and 
[7 = 0. The spin-spin correlation functions at half-filling 
in second-order perturbation theory are given by 



{SiSjj)^Zl-i = 12 ( f )' (j)\f\s^,sI_,\F) (2) 
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where \F) is the ground state of the noninteracting chain. 
This means that for very strong coupling the correlation 
function (S/S//)'-^^ is effectively given by (s2s|^_]^) be- 
tween the first and the last sites of the tight-binding chain 
in-between the two strongly bound Kondo singlets. Such 



FIG. 3. (Color online) DMRG data for (S/S/i}^ on a double- 
logarithmic scale for U — and the impurities at the ends of 
the chain, including a comparison with perturbation theory 
(triangles) for J — —60. 
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FIG. 4. (Color online) Double-logarithmic central differences 
a{R) for U — and various values of J including the results 
of the perturbation theory for |J| ^ 1. 



a behavior can indeed be observed in Fig. [3l where the 
DMRG data for the correlation functions are plotted on 
a double-logarithmic scale and are also compared to the 
perturbation theory for large values of \J\. There is very 
good agreement for J = —60 with relative deviations of 
less than 0.5 %. 

The results from Fig. [3] suggest a power-law behavior 
for large couplings as given by Eq. ([1]) . In order to extract 
the exponent and to validate the existence of a power 
law from our DMRG calculations, we have computed the 
double-logarithmic central differences 

l0g((SjSjj)fl+2) - l0g((S/Sjj)fl-2) 

"^""^^ log(i?4-2)-log(i?-2) 

of the correlation functions for various couplings. They 

will directly give the exponent a if a power law 

is assumed at a certain distance. This analysis is more 



4 




10 100 
Inter-impurity distance R 



2 
1.5 

1 

0.5 



g 1.5 



1 

0.5 



10 



20 



30 



40 



50 



60 



70 



80 



U = 1 



Egg' 



igBBBi 



IEEJ33HEBE 



03 



--X— 

.0... 
□ 



J = — oo 
J = -60.0 
J = -8.0 
J= -2.0 
./ = -0.5 



10 20 30 40 50 60 
Inter-impurity distance R 



70 



80 



FIG. 5. (Color online) DMRG data for (S/S/i)^ on a double- 
logarithmic scale for several values of J at (7 = 1 as well as 
U — 4 and the impurities attached to the ends. 



FIG. 6. (Color online) Double-logarithmic central differences 
a{R) for 17 = 1 as well as 17 = 4 and various values of J 
(J = —ex denoting the perturbative results). The impurities 
are attached to the ends of the chain. 



accurate for the evaluation of the distance dependence 
a{R) than fitting a power law to the data of Fig. [3] In 
Fig. |4]the central differences a{R) are shown for C7 = 
and several values of J. For strong coupling and large 
distances, a fast convergence a — ?■ 2 can be clearly ob- 
served. The asymptotic behavior is also in agreement 
with our perturbative results whose derivative according 
to Eq. (IH) is also shown in Fig. [H For small | J|, the 
accessible distances are too small to find a behav- 
ior. However, the results are in agreement with such a 
behavior for i? — > oo. 

As a curiosity, we remark that there appears to be 
a change in behavior with decreasing |J|. For large 
J|, the asymptotic exponent is consistently approached 
from above, while for smaller values of the coupling, the 
short-distance behavior is weaker than the asymptotic 
We interpret this as a sign of the crossover 
from Kondo screening at large |J|, i.e., outside of the 
scales set by the Kondo correlation length one does not 
have sizable influence by the spins any more, to the 
RKKY dominated regime at small coupling, where the 
correlations for small distances should rather reflect 
the tendency towards antifcrromagnctic order with a 
correspondingly slow decay. 



S. Interacting conduction chain U ^ 

We now switch on the interaction in the conduction 
chain. For |J| — >■ cxd. the system decouples into three 
parts as in the case of a noninteracting chain. Again, 
the impurity-impurity correlation function in this limit 
is proportional to the end-to-end correlation function be- 
tween the two boundary sites of the interacting Hubbard 
chain between the Kondo singlets. For the spin-spin cor- 
relations of the Hubbard chain, we can consult results 
from conformal field theory (CFT). For half- filling, a 



charge gap opens for all C7 > in the Hubbard chain, 
so that the system is a Mott insulator. The spin sector 
on the other hand remains gaplcss. In the continuum 
limit, the Hamiltonian of the Hubbard chain at half- 
filling separates into charge and spinon parts, both per- 
turbed by marginal current-current interactions. From 
the spinon parts logarithmic corrections to the spin-spin 
correlation functions are expected due to the marginal ir- 
relevant perturbation. However, on account of the SU(2) 
symmetry, these logarithmic corrections should have no 
U dependence. From the holon parts we expect expo- 
nentially decaying correction terms e~'^^^^ as two holons 
with mass M are created. These terms are expected to 
decay much faster than the logarithmic correction terms. 
In summary, from conformal field theory for the contin- 
uum limit of the spin-spin correlation functions of the 
Hubbard chain with t/ > 0, we still expect a decay 
with C7-independent logarithmic corrections. However, 
it should be noted that these considerations have to be 
adapted to the special case of the end-to-end spin-spin 
correlations of an open Hubbard chain. 

Hallberg and Egger derived similar results from Lut- 
tinger liquid theory as the spin sector remains gaplcss 
and can therefore still be described by a Luttinger liq- 
uid with g = 0^ With bosonization, they calculated the 
spin-spin correlations for the finite continuous system. 
By evaluating them close to the boundaries, they could 
show that for all values of g, the exponent of the cnd-to- 
cnd spin-spin correlation function asymptotically reaches 
a — ?► 2, at that time in agreement with their DMRG data 
(m = 200). Note, however, that from Fig. [2] we know 
that a too small value for m can dramatically change the 
observed long-range behavior. 

Let us therefore re-examine the long-distance behav- 
ior based on DMRG calculations with m large enough 
to have a reasonably converged behavior at long dis- 
tances. We start by analyzing our numerical data for 
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a pure power-law decay given by Eq. ([T]). Figure [5] shows 
the spin-spin correlation function (S/S//)^ on a double- 
logarithmic scale. Note that this representation suggests 
a very clear power law. The corresponding central dif- 
ferences extracted via Eq. ([4]), which give the exponent 
a, are depicted in Fig. [S] for [/ = 1 and [/ = 4. The 
proportionality from Eq. Q between (S/S//)'^) and the 
correlations (si^L-i) the conduction band has proven 
to be valid for nonzero values of U as well. By calculat- 
ing {s2S^j^_i) for the corresponding open Hubbard chain 
of length L — 2 via the DMRG method the asymptotic be- 
havior of the central differences can be determined by the 
same perturbative approach. These data are labeled as 
J = — oo in Fig. [6l Overall, the DMRG data for nonzero 
U suggest that the exponent a{R) approaches an asymp- 
totic behavior for large couplings and distances, with the 
same difference between strong and weak coupling J on 
how the asymptote is approached. The main result of our 
DMRG calculations, however, is the observation that the 
exponent assumes a value smaller than two. This result 
is, however, in clear contradiction to the expectations 
from CFT and bosonization. 

In order to resolve this contradiction, we now take into 
account the logarithmic correction appearing in the 2kp 
contribution to the correlation function, which is pre- 
dicted by CFT as 

(S,S,,) oc i?-" [log(ii')]"^ . 

From this an expression for the exponent ai of the log- 
arithmic correction is derived via differentiation. One 
obtains 



ai 



dlog(SjSjj) 
d log R 



a log(i?), 



where we now use the exact result a = 2. The expo- 
nent ai as a function of the distance R is depicted in 
Fig. [71 As expected, ai approaches zero rather rapidly 
for increasing distances and U = 0. In the limit J7 — > oo 
the Heisenberg model gives ai — !> 3/2, which fits into the 
picture of an increasing exponent with U. For t/ > 0, 
ai seems to converge to a finite value for increasing i?, 
cf. Fig. [21 Note, however, that even for the largest sys- 
tems studied here, one still finds an increasing value for 
«!. Furthermore, for fixed J7 > 0, ai seems to approach 
the same asymptotic value irrespective of J (not shown 
here). In this sense, our DMRG data are compatible with 
an exponent a = 2 modified by a [/-dependent logarith- 
mic correction. However, the dependence on U is not 
expected from conformal field theory. This contradiction 
will be discussed further in Sec. HVl 



B. Impurities attached to the bulk 

Next, we study the spatial dependence of (S/S//) in 
the bulk limit, i.e., R ^ L/2. For a noninter acting con- 
duction chain in the large-J limit, the system decom- 
poses into five parts: a left and right conduction chain. 
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FIG. 7. (Color online) Exponent Qi of the logarithmic correc- 
tion as a function of the distance R for J = —oo and several 
values of U. 



the two singlets, and an effective chain between them. 
However, the spin-spin correlation functions between the 
two attached spins in this limit arc again just propor- 
tional to the end-to-end spin-spin correlations between 
the two boundary sites of the chain in-between. The left 
and right chains show no dependence on the correlation 
functions in that limit. This was already argued by Hall- 
berg and Eggeri^ Gonsequently, the same expectations 
from conformal field theory and bosonization are valid as 
in the case of the impurities coupled to ends of the chain. 

Due to the oscillatory behavior of the correlation func- 
tion, it is justified to merely consider odd inter-impurity 
distances R in our DMRG calculations. Additionally, we 
choose even system sizes L as this ensures a nondegener- 
ate ground state. For i? = 3, 7, 11, . . . the impurities are 
attached symmetrically around the center of the chain 
whereas for R ~ 5,9,13,... they are shifted out of this 
symmetric setup by one site. This provides us with the 
L-independent values of (S/S//)fl (cf. Appendix [B]) . We 
worked with chains of length L — 120 to make sure that 
the impurities are sufficiently far away from the ends to 
minimize boundary effects. 

Our DMRG results for the correlation functions 
between two impurities attached to the bulk are shown 
on a double-logarithmic scale in Fig. [8l The qualitative 
behavior of (S/S//)^ is very similar to that observed for 
the impurities attached to the ends, cf. Figs. [3l and [5l 
The upper panel of Fig. [9l contains an analysis by means 
of double-logarithmic differences for two impurities at- 
tached to the bulk of a noninteracting conduction chain. 
As for the system with the spins at the boundaries, 
both the large-J and large- i? behavior obey a power law 
with an exponent of a 2 for U = 0. Moreover, the 
results are still consistent with perturbation theory. In 
the presence of correlations U ^ in the bulk limit, we 
still expect the same decomposition of the system into 
five parts. Therefore, the asymptotic behavior for [/ = 4 
in the lower panel of Fig. [SI also turns out to be the 
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FIG. 9. (Color online) Double-logarithmic central differences 
a{R) for U = and ?7 = 4 in the bulk limit for various 
couplings J {J = —CO denoting the perturbative calculation). 



same as for two impurities coupled to the ends in Fig. [6l 
Consequently, if we assume a pure power law the bulk ex- 
ponent assumes a value which is smaller than two as well. 



In Fig. [3] it also has to be noted that the asymp- 
totics are approached faster in the case of stronger 
correlationsi^ This can be made plausible because the 
RKKY interaction, which is mediated via the spin polar- 
ization of the conduction band electrons by the magnetic 
moments of the impurities, is suppressed in the presence 
of stronger correlations. Thus, the crossover from RKKY 
to Kondo behavior occurs for smaller distances and cou- 
plings. 



IV. SUMMARY AND DISCUSSION 

Let us summarize our results for the long-distance be- 
havior of the correlation function between the two impu- 
rities. We have numerically shown that the long-distance 
behavior for the two impurities in the bulk is equivalent 
to the long-distance behavior for two impurities coupled 
to the ends of the chain. This allows us to focus on the 
setup of one impurity coupled to either end of the chain - 
which is a computationally less demanding system. 

It turns out that for a fixed value of U the exponents 
of the power-law decay will converge to the same value 
for different J < 0. The larger |J|, the better the con- 
vergence, which can be seen in Figs. S] and [HI This has 
already been pointed out by Hallberg and Egger.— The 
last open question is whether the exponent depends on 
the interaction U in the Hubbard chain. In order to an- 
swer this question, we can focus on the large- J limit. We 
have shown that our perturbation calculation agrees very 
well with our numerical data. For U ~ 0, wc find a value 
of a = 2 for the exponent which is in agreement with the 
result for the correlations of the noninteracting chain. 
For J7 > 0, we also find seemingly converged exponents, 
but with a < 2. This is in contradiction to expectations 
from CFT and bosonization. In order to resolve this con- 
tradiction, we have analyzed the data for a decay 
with logarithmic corrections. Taking into account the 
logarithmic correction appearing in the 2kp contribution 
to the correlation function, we can show that our results 
are in agreement with a converged exponent a = 2, which 
would be consistent with CFT/bosonization. However, 
we obtain a t/-dependent exponent ai = ai(U) for the 
logarithmic correction, which again stands in contradic- 
tion to the results from CFT. 

The interpretation of this result is difficult. One has 
to remember that the CFT predictions, in particular the 
logarithmic corrections, are based on the continuum limit 
and do not include boundary effects. For the end-to- 
end spin correlations a J7-dependent logarithmic correc- 
tion may thus develop for finite systems if, due to the 
boundaries, a coupling between spin and charge degrees 
of freedom occurs. As one possible consequence, the U- 
dependent exponents of the logarithmic corrections may 
survive in the thermodynamic limit. Since the correc- 
tions are of logarithmic nature, we can infer that one 
naturally needs to access system sizes that are several or- 
ders larger. To our knowledge, such investigations have 
not been performed yet, but would be interesting in view 
of the results presented here. Alternatively, one may 
study a modified Hamiltonian in which the coupling to 
the marginal perturbation is reduced in order to suppress 
the logarithmic corrections. This can be achieved, e.g., 
by adding an additional nearest-neighbor Coulomb inter- 
action V = — C//[2 cos(2fci?)] in the extended Hubbard 
modeli^ 

Furthermore, our observations urge us to be extremely 
careful in extracting exponents from numerical data. The 
double-logarithmic presentation of the correlation func- 
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tion suggests a clean power law and the extracted ex- 
ponents look converged for the system sizes studied and 
suggest a value oi a < 2 ior U > 0. However, we know 
that for much larger systems the true asymptotic value 
of a = 2 with logarithmic corrections will eventually be 
reached. A similar behavior has also been observed for 
the boundary exponent of the spectral function in open 
Hubbard chains. Here, for certain interactions in the 
studied systems of at most 500 sites, the power-law be- 
havior with respect to the expected exponents did not 
occur at all.—"— In the light of these findings, it is clear 
why our studies based on very accurate DMRG simu- 
lations lead to these strong deviations to the observa- 
tion of the seemingly proper, J7-independent exponent 
q; = 2 by Hallberg and Egger. As their calculations 
were based on smaller distances with significantly less 
states kept, we can understand this discrepancy from 
the extreme sensitivity of the power-law exponents on 
the numerical accuracy. In order to fully resolve the 
long-distance behavior of correlations in the two- impurity 
Kondo- Hubbard chain, one would have to analyze system 
sizes that are much larger rendering further calculations 
using the DMRG impossible. In particular, the intrin- 
sic exponential decay of long-range correlations in the 
DMRG makes these calculations of end-to-end correla- 
tion functions very expensive. Possible other algorithms 
to overcome this problem are, for example, MERA^"— 
or other tensor networks, which are known to be able to 
recover the right form of the correlations. 



a single impurity spin can be stated as 
20 / f \ ^ 

(Al) 

where is the ground state of a nonintcracting 

chain of length L — I and V'i^i are fermionic field op- 
erators at site i. If one wants to introduce a second im- 
purity at the other end of the chain, there is a slight mod- 
ification as the first and the last conduction band sites 
will be quenched out. Thus, one uses \Fl-2) instead of 
and also obtains an expression for {Sjj)^^\ which 
is equivalent to Eq. (|Aip with i = L — 1. Combining these 
two expressions yields the impurity-impurity correlation 
function 



(A2) 



in the strong-coupling limit as the two singlets are 
screened from each other by the free chain of length L — 2 
in-between. Here sf = {ipj -f — ipj j^V'i 4.)/2 denote the 
spin operators. The correlations (sls^.i) can be calcu- 
lated analytically for a noninteracting chain at half-filling 
and for S^^^ = 0. Inserting the explicit representation of 
the field operators 



(t) 



L 



{fc„} 
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where fc„ — ni: / [L + 1) and n = 1, . . . , L, gives 



L/2 



(^Lk24-ll^L) 



1)2 5Z ^ 



1=1 q=L/2+l 



M 

sm I nl X 

L + 1 



27r \ fL-1 \ f 2tt 
X sm I — -/ I sin | — ^T^q sin -q 



L+1 



L + 1 



L + 1 



Now Eq. (I3|) directly follows upon exploiting this result, 
Eq. (jA2p and the rotational invariance of the system. 



Appendix B: Bulk setup 



Appendix A: Perturbation theory 

Let us first consider the case of a single impurity cou- 
pled to the first site of a noninteracting chain. In the 
strong-coupling limit |J| ^ 1, there is hardly any hop- 
ping between the sites 1 and 2 due to the rigidly bound 
singlet between the impurity and the spin at the first con- 
duction band site. So, the corresponding hopping term 
represents the perturbation. The calculations in Ref. [lol 
reveal that the expectation value to second order in t of 



As described in Sec. IIIIBl we regard odd distances R 
and even chains of length L to ensure a nondegenerate 
ground state. If one attaches the impurities symmet- 
rically around the center of the chain, the correlation 
function will have nearly the same value for {R + l)/2 
even irrespective of L as long as R ^ L/2 is fulfilled, 
while it will have a strong size dependence for other odd 
values of R. This is shown for J7 = 4 and J = — 2 in 
Fig. [TOl By means of a finite-size scaling with respect 
to 1/L the extrapolated values for the other distances 
i? = 5,9,13,... will obey the same behavior for L — )■ cxo 
as the L-independent i?'s. 
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FIG. 10. (Color online) Finite-size effects in the bulk limit 
[R <Si L /2) for impurities attached symmetrically around the 
center of the chain for U — 4 and J = —2. 
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